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Abstract 

We analyze the charm quark electric dipole moment (EDM) in the minimal supersymmetric ex- 
tension of the standard model (MSSM), including important two-loop gluino contributions. Con- 
sidering the experimental constraint on the neutron EDM, the theoretical prediction for the charm 
quark EDM can reach about lO'^'^e • cm. If taking into account the mixing between the scharm 
and stop quarks in the effective supersymmetry scenario, the charm quark EDM can be enhanced 
to 10~^^e • cm. Direct production of the CP-odd ^Pi state he in e^e~ annihilation via the 
CP-violating process at the BES-III and CLEO-C is analyzed. 

PACS numbers: ll.30.Er, 12.60.Jv,14.80.Cp 
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I. INTRODUCTION 

Up to now, CP violation is only found in the and 5-systemj^, which can be well 

explained within the standard model (SM) of electroweak interaction. As is well known, 
electric dipole moment (EDM) of an elementary particle is a clear signature of CP violation. 
However, EDM of a fermion does not appear up to two-loop order, and the three-loop 
contributions partially chancel among each others in the SM. Therefore, observation of 
sizable EDM of an elementary fermion would definitely be a signal of existence of new 
physics beyond the SM with extra CP phases. 

Supersymmetry (SUSY) is now believed to be the most attractive scenario for new 
physics. Beside the Cabibbo-Kobayashi-Maskawa (CKM) mechanism, the soft breaking 
terms provide a new source of CP and flavor violation in the minimal supersymmetric 
extension of the standard model (MSSM). Those CP violating phases can affect some im- 
portant observables in the mixing of Higgs bosons [3], the lepton and neutron electric dipole 
moments (EDMs) , lepton polarization asymmetries in the semi-leptonic decays 16] , the 
production of P-wave charmonium and bottomonium and CP violation in rare i?-decays 
and B^B^ mixing At present, the strictest constraints on those CP violation phases orig- 
inate from the lepton and neutron EDMs. Since the lepton and neutron EDMs have not 
been measured so far, there are several suggestions to keep the theoretical estimate of the 
neutron and electron EDMs below the experimental limits, they are (i) choosing small CP 
phases < (9(10~^)[^, (ii) finding appropriate parameter domain where various contributions 
cancel with each other [lO|, or (iii) making the first two generations of scalar fermions heavy 
enough (heavier than 20 TeV) while keepin g th e soft masses of the third generation below 
TeV to keep the Higgs boson naturally light 

In this work we calculate the charm quark EDM in the CP violating MSSM. As pointed 
out in 4], the quark chromoelectric dipole moment (CEDM) makes important contributions 
to the quark EDM at a low energy scale. In addition, there is a parameter space where 
some two-loop gluino diagrams provide important contributions to the quark EDM. In our 
analysis, we include all the contributions from those pieces. We first work in the simplest 
model, neglecting the new possible flavor violation sources except the CKM mechanism, 
and avoiding any assumptions about unification of the soft breaking parameters. We find 
that the charm EDM can reach about 10~^''e ■ cm. Then we also calculate the charm quark 
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FIG. 1: The one- loop self energy diagrams which lead to quark EDM and CEDM in the MSSM, the 
corresponding triangle diagrams are obtained by attaching a photon or gluon line to the self-energy 
blob in all possible ways. 



EDM in the effective supersymmetry scenario ilfij where the scenario (iii) mentioned above 
is accommodated. Since the first two-generation sfermions are very heavy, the charm quark 
EDM is greatly suppressed in this scenario. However, while considering large mixing between 
the second and the third generation sfermions we find that the charm quark EDM can be 
greatly enhanced to as large as about 10"^^e ■ cm. The reason is that, since the effective 
operator which determines the charm quark EDM induces a chiral flip, the large mixing 
between the scharm and stop can enhance the charm quark EDM due to the top quark 
mass. 

The analysis in Refs. 
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18| shows that the EDM of heavy quark plays an important 
role in the direct production of singlet P-wave quarkonia, thus with upgraded BEPcjlfl| 
aud CLEO-C programsW. it might be possible to produce the 'P, charmonium state, ft., 
directly. Conversely, if no he meson is observed, an upper bound on the charm quark EDM 
would be set. 

The paper is organized as follows: In the next section, we give the analytic expressions 
of the charm quark EDM and CEDM in CP violating MSSM. In Section III we present the 
numerical results. Section IV is devoted to discuss he production which is closely associated 
with the magnitude of the charm quark EDM, at e^e~ colliders. Finally we conclude in 
section V. 



II. THE CHARM QUARK EDM 



In the effective Lagrangian, the charm quark EDM is defined through a dimension-five 
operator 

I^Euu = -ldJ:a^^^,cF^^ (1) 

with F^^ being the electromagnetic field strength. In a theoretical framework with CP 
violation, the corresponding loop diagrams induce the charm EDMs. Since quarks also take 
part in the strong interaction, the chromoelectric dipole moment (CEDM) operator is 

where T° (a = 1, ■ ■ ■ , 8) denote the generators of the strong SU{?>) gauge group and G^^, 
is the chromo-electromagnetic field strength. This operator contributes to at low energy 
scale. In principle, the dimension-six Weinberg operator and two-loop Bar-Zee diagrams 
all provide contributions to d^. Nevertheless, the contributions from the Weinberg operator 
and Bar-Zee diagrams to d^ can be ignored safely as we consider the constraint on the 
neutron EDM. The best way to describe certain loop-induced contributions is the effective 
theory approach, where the heavy particles are integrated out at the matching scale and 
the effective theory includes a full set oi CP violating operators. In this work, we restrict 
ourselves to the following operators that are relevant to the c-quark EDM 

C = -lrf7c^A*-^5cF - -rf3cTV'^'^75cG" . (2) 

The one-loop supersymmetric contributions to the Wilson coefficients in Eq. Q originate 
from three types of graphs: gluino-squark, chargino-squark, and neutralino-squark loops 
(FIGH)). Details about the descriptions are presented in ref.|l3|. The contributions of the 
one-loop gluino-squark diagrams are 
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Here a — e^/ (47r), 9^ denotes the phase of m^, and Z- {q — u, - • ■ ,b) are the mixing matrices 
of squarks, i.e. Zjm^Z, = diag{m^^, m^^) where 



m2+m; + m2(i-Q^s2) cos 2/3 m^(A*-//i?J 

mfA -/x*i?J m2_ _ +m2 + m2(l-g s2)cos2/3 

with = 2/3 (—1/3), i?^ = tan/3(l/ tan/?) for g = -u (d). tan/3 = ^ is the ratio between 

d 

the up- and down-type Higgs vacua, and is the Weinberg angle. The shortened notations 
— sin^^, = cos^^ are adopted. The loop functions are 



B{r) = [2(r-l)2]-^[l+r+2rlnr/(r-l)], C(r) = [6(r-l)2]-^[10r-26-(2r-18) lnr/(r-l)]. 
In a similar way, the one loop neutralino-squark contributions can be written as 
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with 
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X (^g)i,i(^iv)4,fe ; 
m c 



. (6) 

Here OL^—g\l (47r), = sin /3, = cos /3, and is the mixing matrix of neutralinos. 
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Finally, the chargino-squark contributions are formulated as 




FIG. 2: The two-loop self energy diagrams which lead to quark {q = u (d), Q = d (u)) EDM and 
CEDM in the MSSM, the corresponding triangle diagrams arc obtained by attaching an external 
photon or gluon line to the self-energy diagrams in all possible ways. 
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where V denotes the CKM matrix, and the loop function is 

A{r) = 2(l-r)~2[3-r + 21nr/(l - r)]. 



The couphngs are defined as 



Tfl 
Tfl 

where Z± are the right- and left-handed mixing matrices of the charginos. Noting that 
one-loop chargino contributions to the quark EDMs and CEDMs are proportional to a 
suppression factor ^ which exists in all the couplings (g = m, d). 

The two-loop gluino contributions originate from the two-loop self energy diagrams (FIG. 
12)). The corresponding triangle diagrams are obtained by attaching an external gluon or 
photon line to the self-energy diagrams in all possible ways. In those two-loop diagrams, 
there are no new suppression factors except the loop integration factor in the corresponding 
amplitude, and the resulting Wilson coefficients of those dipole moment operators do not 
involve ultra-violet divergence. Furthermore, there is a parameter space where the two-loop 
results are comparable with those one-loop contributions because the dependence of the 
two- loop results on the relevant CP phases differs from that the one- loop results on the 
corresponding CP phases. The contributions of the two-loop gluino-gluino diagrams to the 
quark EDMs and CEDMs are formulated as 

8e ea^lm, 

d^ 



5(2) 9(47r)2m: 



-FAx ,x- )Imf(Z-), .(Z^)^^e~*^3 

9(2) 9(47r)2m2 ^ ® V ''^^'•'^ s'l'^ / ' ^ ' 

where Xi = mf/m^, and the function Fi{xo, xi, X2, xs, X4) {i = 1, 3) are defined in Ref . [l^ . 
The contributions of the two loop neutralino-gluino to the quark EDM and CEDM are 
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given by 



q s 



< |m-|F4(x ,x- ,x-,x „,x- ) Im( (A^ )..(Z^)., 
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As for the two-loop gluino-chargino contributions to the c-quark EDM and CEDM, we 
have 



Aeaa, 



-VlVn \\m-\FG(x^,x- ,x-,x +,x-)\lm((A^)^.(Z^)., 
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Note that the last terms of the (P, and d^, are not proportional to the suppression factor 

X±(2) xf(2) 

m^/m^. This implies that the two- loop gluino-chargino diagrams may be the dominant part 
of the chargino contributions to the quark EDMs and CEDMs. 

In order to account for resummation of the logarithmic corrections, we should evolve 
the coefficients in the quark EDM and CEDM operators at the matching scale ii down to 
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the charm quark mass scale ~ 1.2 GeV jlj| using the renormahzation group equations 
(RGBs) 



(12) 



where r]^ ^ 1.53 and rj^ ~ 3.4. At the low energy scale, we need to include the contributions 

from the quark CEDMs to the quark EDMs, when we evaluate the numerical value of the 

I — I 

charm quark EDM. This is realized by a naive dimensional analysis as 



III. NUMERICAL RESULT 

A. In case of no flavor mixing 

In this subsection we present our numerical results where no flavor mixing is introduced 
in the squark sector. At present, the upper bound on the neutron EDM is (i„ < 1.1 x 
10~^^ e ■ cm. In order to make consistent with this constraint, we take the squark masses 
of the first generation as = = = 20 TeV, the CP phase of the /i parameter 
9^ = 0. Without losing generality, we choose the absolute values of the /i parameter and 
soft trilinear quark couplings as = A^^ = 300 GeV with q = u, d, ■ ■ ■ , t, and the soft 
SU{2) X U{1) gaugino masses as |mj = \m^\ = 600 GeV. For the soft masses of the third 
generation squarks, we set m, = m. = m, = 500 GeV. 

Taking tan/? = 2,arg{A^) = 0, the SU{3) gaugino mass \m.\ = 400 GeV and the soft 
masses of the second generation squarks as m, = m. = m, = 300 GeV, we plot the 
c-quark EDM d^ versus the CP phase of m. in FIG. 0] With our choice of the parameter 
space, the c-quark EDM d^ can reach 0.5 x 10^^'^ e ■ cm while the neutron EDM d^^ is below 
the experimental limit. Numerically, we find that the dominant contributions to d^ originate 
from the gluino-squark sector (including the one- and two- loop contributions). The sum of 
the char gino- squark one-loop and chargino-gluino- squark two-loop contributions to d^ is less 
than 10~^^ e ■ cm. Through the CKM mechanism, the contributions to d^ from the second 
squarks (including squark- char gino one loop and squark- chargino-gluino two-loop diagrams) 




(13) 
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FIG. 3: The c-quark EDM varies with the CP phase of m. at tan (3 = 2, where the dashed Une is 
the result at the one-loop order, and the solid line is the result at the two-loop order. The other 
parameters are taken as shown in text. 

is below 10~^^ e ■ cm. This fact can help us understanding why satisfies the experimental 
bound, while can reach a relatively large value. The choice of the parameter space in FIG. 
E]is the same as that in FIG.IHlexcept for tan/? = 20. Comparing with FIG.|2l the maximum 
of d^ in FIG. 01 is about 0.9 x 10~^° e ■ cm. Since the dependence of gluino-squark two-loop 
contributions to the quark EDMs on the CP phases is the same as that of gluino-squark 
one- loop contributions to the quark EDMs on the CP phases, the gluino-squark two-loop 
contributions to d^ is below 10% with our choice of the parameter space. 
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FIG. 4: The c-quark EDM varies with the CP phase of m. at tan j3 = 20, where the dashed hne 
is the result at the one-loop order, and the solid line is the result at the two-loop order. The other 
parameters are taken as in text. 

Now, we investigate how varies with the CP phases of the soft squark Yukawa cou- 
phngs. Taking tan (3 = 2, arg{m-) = 0, we plot versus the CP phases arg{AJ = arg{A J 
in FIG. El where the other parameters are the same as in FIG. El It is clear that the maximum 
of the absolute value of d^ is about 0.9 x 10~^° e-cm. With the same choice of the parameter 
space as for FIG.Elexcept tan/3 = 20, d^ can reach 5x 10^° e-cm. Nevertheless, the chargino 
sector which induces a contribution to d^ is above the experimental upper limit. In order 
to suppress the chargino contributions to , we choose m. = m. = m . = 500 GeV. 




11 



-100 I I 

0.0 0.5 1.0 1.5 2.0 

arg{A^)/'n: = arg{AJ/Tr 

FIG. 5: The c-quark EDM varies with the CP phases of soft squark Yukawa couphngs arg{A^) = 
arg{A^) at tan/3 = 2, where the dashed hne is the result at one-loop order, and the solid line is 
the result at two-loop order. The other parameters are taken as in text. 

Correspondingly, we set the SU{3) gaugino mass as |m-| = 600 GeV. Within this scenario, 
the two-loop gluino-squark contributions are much less than that of one-loop gluino-squark 
contributions (FIG. Ell- 
in a similar way, we can investigate how varies with the CP phases of the SU (2) x U{1) 
gaugino masses. However, the contributions from the phases of the SU{2) x U{1) gaugino 
masses to are below 10^^^ e ■ cm for our choice of the parameter space. 
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argiA^)/^ = arg{AJ/ir 

FIG. 6: The c-quark EDM varies with the CP phases of soft squark Yukawa couphngs arg{AJ = 

arg{A^) at tan/3 = 20, where the dashed hne is the result at one- loop order, and the solid line is 
the result at two-loop order. The other parameters are taken as in text. 

B. Effective SUSY 

In this subsection, we consider the case of the effective SUSY scenario including mixing 
between the second and the third generation squark. We fix the soft squark masses of the 
first two generations as heavy as 20TeV. Wc denote the different left-handed and right- 
handed squark masses as rni^ff^^ ^ ^ respectively and parameterize the the general form of 
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l.e-20 



mi: =mR =20 TeV, mj: =mR = 300 GeV 
0L=9iR=6, A=600i, my2=26o GeV 




-0.2 0.0 0.2 
(Or) M 



FIG. 7: dc as a function of 9i = 9ji with tan/3 = 5, ^ = 150Gey, m 
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m 
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20TeV, 



m 



= = SOOGeV, mi/2 = 200GeV, A = 600iGeV , m-g = 2hf)GeV and ./.^ = 0r = 0. 



the soft breaking mass matrix as 



(14) 



with 



v 



cosOl 
sin eLe-'"t'^/'^ 



sin Olc''^^/'^ 



( 2 



, and m| 



2 

^2 



V 



m 



(15) 



is / 



and the form for the right-handed sector is similar. The free parameters are now tan/?, mi^, 
^R^y ^L,R, 'pL,R, and A = \A\e^'^^ for the squark sector; = \^\e'"^^' and mi/2 for the chargino 
and neutrahno sector; and the gluino mass m^. Assuming the gaugino masses are universal 
at the GUT scale, the phases of the gaugino masses can be rotated away. Compared with 
the regular cases, the effective SUSY scenario has four additional parameters: 6l,r, 4>l,r- 

By our calculation, we are convinced that the dominant contribution to dc comes from 
exchange of gluino. Therefore the parameters related to the down squark sector and the 
chargino and neutralino sector do not affect the numerical result significantly. This fact helps 
us to reduce the parameter space. We always take the same value of mi_, 9l,r and (j)L,R for 
the up and down squark sector. We set tan (3 = 5 and /i = 150GeV in our calculations. 

By the above discussions we are left with only 8 free parameters: mi^{=mj^^), 6l,r, 4>l,r, 
and A = l^le*"^^. We first show the effects of the mixing angle 6l,r in Fig. [7| We notice 
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FIG. 8: dc as a. function oi 6^ (long dashed line), Or (short dashed line), 6l = —Or (solid line) 
and Ol = On (dotted line) for 0^ = Or = 7r/4. Other parameters are the same as those in Fig. EJ 
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FIG. 9: dc (left panel) and the branching ratio of t — > 07 (right panel) as functions of |^| for 
maximal mixing angle = = 7r/4 and (j)L,R = 0. Other parameters are the same as those in 

Fig. m 



that if the mixing angle is zero, dc is nearly zero, which is a direct result of extremely heavy 
charm squark of about 20TeV. For maximal mixing, 6l,r = vr/4, we get the maximal dc- 

The contribution to dc from the phases of (f)L,R is plotted in Fig. |H1 for maximal mixing 
angle 9l = Or = 7r/4. It is interesting to notice that when 9^ = Or (dotted line) the 
contributions from 0^ and (f)R cancel each other. 

The dependence of dc on the magnitude of A is plotted in Fig. 121 (on the left panel). For 



15 



7 



'Pi 



Z 

(a) 



'Pi 




FIG. 10: The Feynman diagrams for production of the ^Pi charmonium resonance in e~^e~ scat- 
tering from (a) the CP-conserving and (b) CP-violating modes. The dark blob stands for insertion 
of the effective operator 

large |y4|, not only the CP violation is enhanced, but also the effects of the lighter stop mass 
is suppressed. Therefore dc increases rapidly with increasing value of \A\. The maximal 
value of dc can reach about ~ lO^^^e ■ cm. 

The numerical results show that in the effective SUSY scenario, where all the bounds set 
by the low energy experiments are satisfied, the EDM of charm quark can be enhanced by 
an order of magnitude due to the large mixing between the second and third generations of 
up squark. 

Finally, we show the flavor changing effects of large mixing between scharm and stop. 
The branching ratios of t ^ 07 and t ^ eg induced by the mixing are plotted on the right 
panel of Fig. El which are below ~ 10~^ and ~ 10~^ respectively. We have adopted the 
LoopTools package 2l| to finish the loop integrations. In fact, the present experiments do 
not constrain the mixing between scharm and stop 
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IV. hr PRODUCTION IN e+e COLLIDER 



Recently the E835 Collaboration at Fermilab |23| and the CLEO Collaboration 24 1 have 
announced the observation of the ^Pi CP-odd state he- However, the direct production of 
the 'Pi charmonium resonance at the e"'"e~ annihilation plays an important role in probing 
the charm quark EDM, since the coupling of photon to the 'Pi charmonium, as shown in 
Fig. Uni (b), is identical to the effective operator in Eq. ((T)). The coupling of the current 
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density [2 

J^(cc|^Pi) = c(a;)(9^75c(x) (16) 

to a photon results in the effective operator (0). We can also check; the quantum numbers 
of the process in Fig. ^l(b). Transition from the initial ^Si {J^^ = 1 ) state to the final 
^Pi {J^^ = 1"*"") state requires the photon-charm-charm vertex, denoted as dark blob in 
Fig. ^1 (b), to be — or + under P and C transformations. We find that the operator 
just satisfies the requirement. Under P and C transformations we have 

The operator (1) obviously possesses right quantum numbers. 

In the SM, the h^. (1^-Pi) meson can be produced via 7Z and ZZ box diagrams, as shown 
in Fig. El (a), which are CP-conserving processes. As shown in Ref. jl?!. \\\. the effective 
Hamiltonian is given as 

HsM = ^^G',.mem,fBJ^(cc|^Pi) ■ J^{e+e-\'P{) , (19) 

where the current is defined in Eq. and 03 is a standard loop integral. In the minimal 
SUSY, there is another diagram where Z is replaced by a CP-odd Higgs and generally 
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it is considered to be smaller than the SM contribution 

The effective Hamiltonian responsible for the CP-violating process of he production via 
the charm quark EDM, shown in Fig. El(b), is given as 

HsusY = ^-J,{cc\'Pi) ■ J^(e+e-|=^5i) , (20) 

where J^(e"'"e~ 1^5*1) = e'y^e. Comparing the two contributions we obtain the production 
amplitude of hf. via the CP-violating process and it can overtake the CP-conserving one if 
the EDM of the charm quark exceeds the critical value [l^ 



log — ~ IQ-^^cm . (21) 



12V27r2 M| ^ m. 

The numerical results presented in last section show that the EDM of the charm quark in 
the effective SUSY is well above the critical value 10~2^e ■ cm. Therefore, direct production 
of he at e^e~ annihilation is expected to be closely related to the EDM of the charm quark 
or even determined by it. 
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Assuming that the center-of-mass energy of BEPC or CLEO-C is set at the vicinity of 



MhX^ 3524Me\^) 2J|, we get the cross section for production as 
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a(e+e" ^ K) = 27 



dc 
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^5(0) 



cr e e 



Si. 



27 



Mhjc 



(Tee 



Si, 



(22) 



where Rp and Rs are the wave functions of he and ^Si states, respectively. The smallness 
of the EDM of charm quark predicts a very small he production rate, which is about 10 
orders of magnitude smaller than a(e~^e^ Si) for h^ ~ lO^^^e ■ cm. With the upgraded 
BEPcQ and CLEO-C programsQ, a data sample of 10^ ~ 10^^ J/^ will be collected. We 
can estimate the number of hr events asl 



.0' 



Nh, 27 X Nj/^\ ^—\ ~ 



10 



-10 



dc/e 



10 



-20 



cm 



X 



(23) 



which predicts maximally 10^ ~ 10^ h^ events per year in the effective SUSY scenario with 
the upgraded BEPC and CLEO-C luminosity. With the high energy resolution detector 
at BES-IIlj26|, the spectrum of gamma at the inclusive channel he ^i]c can be well 
with high efficiency. The directly produced he may be in the observable range 



determinec 
at BES-III 



V. CONCLUSION 



Considering the neutron EDM constraint, we investigate the EDM of charm quark in 
the CP violating MSSM. Typically, d^ can reach about 10^^° e ■ cm as we include the 
contributions from c-quark CEDM. The mixing between the second and third generation 
squarks in the effective SUSY scenario further enhances d^ to about 10~^^e ■ cm. 

Large EDM of charm quark can induce a direct production of the l^Pi CP-odd state 
charmonium, he, via a CP-violating process at e"'"e~ annihilation. By the luminosity and 
efficiency of the upgraded BEPC-II or CLEO-C (with lO^^J/ip collected) if no he is observed, 
a rigorous constraint on the charm quark EDM would be set, as de < lO^^^e ■ cm. On the 
contrary, any signal of he production will indicate large EDM of the charm quark and 
becomes a clear evidence of new physics beyond the SM. 

The BEPC-II and CLEO-C will provide excellent opportunities to probe many important 
theoretical objects. Among them, the charm quark EDM is a prominent one because it is 



a clear signal of CP violation and an evidence for new physics beyond the SM 
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Our 



numerical results indicate that the direct production rate of he at the vicinity of Mh^ 



2,b2AMeV 



24| may be dominated by the charm quark EDM as long as new physics such as 



the MSSM is involved. Observation of h^. thus produced would clearly indicate new physics 
beyond the SM, conversely, a negative result would set a rigorous constraint to the parameter 
space for the new physics. Therefore we are looking forward to the new data of BEPC-II 
and CLEO-C to help us gaining insight to the new physics beyond the SM. 
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